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Gauge-independence of gluon spin in the nucleon and its evolution 
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In recent papers, we have established the existence of gauge-invariant decomposition of nucleon 
spin, each term of which can be related to known high-energy deep-inelastic-scattering observables. 
A subtlety remains, however, for the intrinsic spin part of gluons at the quantum level. In fact, it 
was sometimes claimed that the evolution of gluon spin depends on the gauge choice and its physical 
interpretation makes sense only in the light-cone gauge. In the present paper, we will demonstrate 
explicitly that the gluon spin operator appearing in our decomposition evolves gauge-independently 
and that it properly reproduces the familiar evolution equation for the 1st moments of polarized 
quark and gluon distributions obtained with the Altarelli-Parisi method, which cannot directly be 
checked by the standard operator expansion method. 



In analyzing the decomposition of the nucleon spin, 
color gauge-invariance is still the cause of intense debate 
[H -[l6|. Although it has long been believed that the 
total angular momentum of the gluon cannot be gauge- 
invariantly decomposed into its intrinsic spin and orbital 
parts [l|, [4[, t his mistaken belief has recently been ques- 
tioned [9[-[l^. In fact, we now have two different gauge- 
invariant decompositions of nucleon spin into the spin 
and orbital angular momenta of quarks and gluons. The 
one is the decomposition proposed by Chen et al. [§| , 
and the other is the one advocated by the present au- 
thor [ll|,[12|. A crucial advantage of the latter decom- 



position is observability of each term of decomposition 
by means of known high- energy deep-inelastic-scattering 
(DIS) measurements [12| . 

The circumstances for the quark part should already 
be clear from the pioneering works by Ji Hi- First, 



the decomposition of the nucleon spin into the quark 
and gluon total angular momentum, Jq and Jg, is man- 
ifestly gauge-invariant, and each term can in principle 
be extracted from generalized-parton-distribution (GPD) 
analyses. Since the intrinsic quark spin (or the longitu- 
dinal quark polarization) AE is clearly gauge-invariant 
and measurable, the quark orbital angular momentum 
(0AM) defined as Lq = Jq — ^ AE should also be gauge- 
invariant. An important fact here is that the quark 0AM 
obtained in that way corresponds to a nucleon matrix el- 
ement of dynamical quark 0AM operator containing the 
full gauge-covariant derivative (3|.|l7|. not the canonical 
0AM or its nontriviai gauge- invariant extension as advo- 
cated by Chen et al. |9j.|10|. 

Although slightly more delicate, the situation for the 
gluon part appears quite analogous. It is widely believed 
that the gluon spin (or the longitudinal gluon polariza- 
tion) Ag is a measurable quantity from the polarized DIS 
experiments. Then, if one defines the gluon 0AM by 
Lg = Jg — Ay, Lg \s clcarly a measurable quantity. A key 
question here is whether Ag is really a gauge-invariant 



quantity or not. 



anticipated), it is a logical conclusion that i-ig 



If it is indeed the case (as is vaguely 

Lg is also 

gauge-invariant. In this sense, a convincing check of the 
gauge-invariance of gluon polarization in the nucleon is 
a fundamentally important homework left in the nucleon 
spin physics. The gauge-invariant glu on spin operator 
was first proposed by Chen et al. and it was later 

confirmed by us 11 [Jl2| . Our investigation goes further. 



We could show in [12[ that the gluon 0AM Lg defined by 
Lg = Jg — A (7, or more precisely the difference between 
the 2nd moment of GPD Hg{x, ^, t) + Eg{x, ^, t) and the 
1st moment of the polarized gluon distribution Ag{x), co- 
incides with the nucleon matrix element of the dynamical 
gluon 0AM containing the potential angular momentum 
not the canonical 0AM [11 



or its 



term explained in [11 [ 
nontriviai gauge- invariant extension advocated by Chen 
et al. [9|,[ig 

In spite of the above-mentioned nice correspondence 
between the quark and gluon sectors, there still remains 
a subtlety in the gluon part. It is no wonder that the 
gluon spin operator given in [l2l | is gauge-invariant at 
least formally, or at the classical level. It was also verified 
in [12] that its nucleon matrix element reduces to the 1st 
moment of Ag{x) in the LC gauge. If this is the case, a 
general thinking based on the gauge pinciple dictates the 
following. Since our gluon spin operator (although not 
necessarily local) is gauge-invariant, its numerical value 
should be the same also in any other gauges than the 
LC gauge. Furthermore, its Q^-evolution should also be 
gauge-independent, as long as the regularization main- 
tains gauge-invariance. 

Unfortunately, life is not so simple. In fact, in an in- 
fluential paper [Si], Hoodbhoy, Ji, and Lu claimed that 
the gluon spin evolves differently in two gauges, i.e. the 
LC gauge and the Feynman gauge (a typical covariant 
gauge). According to them, the calculation in the Feyn- 
man gauge does not reproduce the well-known Altarelli- 
Parisi evolution equation for AT, and Ag [is'l . Undoubt- 
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edly, the difficulty is connected with the widely-known 
observation that there is no gauge-invariant twist-2 lo- 
cal operator corresponding to the gluon spin. A state 
of disorder has been strengthened further by the recent 
claim by Wong et al. [iJl. (See also a similar assertion 
by Cho et al. for the evolution of quark and gluon mo- 
mentum fractions [l^l-) They claim that the gluon spin 
in their decomposition evolves differently from the stan- 
dard Altarelli-Parisi equation. This sounds strange to 
us, because their definition of gluon spin can be thought 
of as a gauge-fixed form of our more general (or formal) 
gauge- invariant expression given in the paper 1^. If 



gauge-invariance is maintained at every stage of manip- 
ulation, the answer should be the same as that obtained 
in the LC gauge. 

Now, the purpose of the present paper is to resolve the 
above-explained puzzle. To this end, we shall explicitly 
calculate the 1-loop anomalous dimension of the quark 
and gluon spin operators appearing in our decomposition 
within the Feynman gauge and show that the answer 
precisely coincides with that obtained in the LC gauge. 

Among the four terms in the decomposition we 
concentrate here on the quark and gluon spin operators 
given by 



M, 
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(1) 

(2) 



The gauge-invariance of M^_^p^^ is self-evident, while 
that of M^_}'^p-^ is ensured by the covariant trans- 
formation property of the physical part of the gluon 
field under a gauge transformation, i.e. ^physi^) ~^ 

We start with writing down a little more explicit form 
of our gluon spin operator : 



- = Va 

g—sptn ^ 



Vb 



with 



= {d'^ A]^) A^a,phys ~ (^^^a) ^a,p/iys i 
"^8 = — [{d^^ ^t) ^a,phys ~ {'^'^ -^a) ■^\,phys _ 

— 9 fabc A^ ( Al Al pf^y^ — Al Al pf^y^ ). 



(3) 

(4) 
(5) 
(6) 



In the LC gauge (^4+ = 0), which falls into the category 
called the physical gauge, only the vertex Va survives 
that simplifies the work process significantly 20]. Note 
however that, in general covariant gauge including the 
Feynman gauge, the gluon spin operator is not given by 
Va alone, but it is a sum of the three pieces Va, Vb and 
Vc- We point out that Vb and Vc have basically the 
same form as the operators Oi ~ — J d^x x Aa and 
O2 = — J d^xfabcA'^Ab X Aa, which, together with the 
naive gluon spin operator 5+ in the LC gauge, were con- 
sidered in the Feynman gauge calculation of the gluon 
spin evolution equation [5|. However, there is a very 



dehcate difference between Va,Vb,Vc in Eqs.Q-dH]) and 
5+, Oi, O2 considered in One should notice that one 
of the gluon fields in our operators Va,Vb, Vc is its phys- 
ical component. (Note that this is an important factor, 
which ensures the gauge-invariance of our gluon spin op- 
erator.) 

The question is then how to introduce this unique fea- 
ture of our gluon spin operator into the Feynman rules 
for evaluating corresponding anomalous dimensions. Re- 
member first that the gluon propagator in general covari- 
ant gauge is given by 



fc2 -I- ie 



with ^ being an arbitrary gauge parameter. The Feyn- 
man gauge corresponds to taking ^ = 1, while the Landau 
gauge to ^ = 0. Since one of the gluon field appearing 
in our gluon spin operator is its physical part, we must 
replace the gluon propagator by 



i5ab 



fc2 



2 



'(fc,A)e''(fc,A), 



(8) 



when one of the endpoint of gluon propagator is obtained 
by the contraction with the physical part of A^ in our 
gluon spin operator. Here, we need a sum of the product 
of gluon polarization vector over two physical polarization 
states (not including the scalar and lon gitu dinal polariza- 



tion). The answer is well-known [21|,[24|. Given below 
is its covariant form, which holds in an arbitrary Lorentz 
frame : 



T'^-EE^e'^(fc,A)£''(fc,A) 



A=l 



n ■ k 



fc^ k" 
{n ■ fc)2 



(9) 



where n being an arbitrary four-vector subject to the 
condition n-e = and n-k ^ 0. For practical calculation, 
it is convenient to take n to be a lightlike four-vector 
with rv^ = 0. In this case, the modified gluon propagator 
reduces to 



k'^ + ie 

iSab 

k^ + ie 



A=l 



k^' n" + k" , , , 
+ n-k (^0) 



which precisely coincides with the gluon propagator in 
the LC gauge. This does not mean that we are working 
in the LC gauge from the very beginning. In fact, if we 
did so, there would be no contributions to the anomalous 
dimensions from the operators Vb and Vc- As pointed 
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out before, it is crucial to use the above propagator only 
when one of the endpoint of the gluon propagator is ob- 
tained by the contraction with the physical part of v4^t in 
our gluon spin operator. In other places, we should use 
the standard gluon propagator, which, for instance in the 
Feynman gauge, is given by 

D'Jik) = j^TTei-a"')- (11) 

The momentum space vertex operators for the gluon 
spin, which takes account of the above subtlety, can be 
expressed by the following formulas supplemented with 
the diagrams shown in Fig IT] : 

- (m^^), (12) 

~ (m^^), (13) 

Vc = 9 fabc ( 9"' 9"' - .9^' 9"' ){P^ +P^) 

+ 9 fabc .9^+ ( 9"' - .9^^' 9^' )iP^ + P^) 

+ 9fabc9-'^i9^'9''' - 9^'9'''){Pt + P^)- (14) 

Here, P^ is a kind of projection operator, which reminds 
us of the fact that we must use the gluon propagator 
_D^^(fc) given by pop . whenever it contains the Lorentz 
index v. 

Keeping in mind somewhat nonstandard Feynman rule 
explained above, we are now ready to calculate the 1-loop 
anomalous dimension of the quark and gluon spin oper- 
ators. The graphs (a),(&),(c) and (d), shown in FiglSJ 
respectively contribute to the 1-loop anomalous dimen- 
sions A7^^\ A7^Q , A7q^ , and ^Jqq- Graphs with ex- 
ternal self-energies are not shown in the figure. Graphs 
which are not symmetric with respect to the vertical lines 
through the operator vertex have to be counted twice. 

For the quark spin operator, we obtain just the ex- 
pected answer : 

a7ig= 0, (16) 

with as — (7^/(4 7r) and Cp = 4/3. Here, the 2nd term 
of A^if comes from the quark field-strength renormal- 
ization in the Feynman gauge. For the gluon spin part, 
we have 

,^(0) Ois ^ as I as 3 

where the 1st term is the contribution of the vertex Va, 
while the 2nd is that of Vb ■ Finally, we find that 

. (0) _ as II ^ as ( 23 \ 

27r 2 27rV6 3V 



with C'a — 3 and riy being the number of quark flavors. 
Here, the 1st, 2nd, and the 3rd terms are respectively the 
contributions from the operators Va , Vb and Vc , whereas 
the 4th term comes from the gluon field-strength renor- 
malization in the Feynman gauge. 

One confirms that the final answer just coincides with 
the lowest-order anomalous dimensions corresponding to 
the longitudinally polarized quark and gluon distribu- 
tions obtained with the Altarelli-Parisi method. We em- 
phasize that, in the standard operator expansion frame- 
work, a direct confirmation of the above answer was done 
only in the LC gauge, because there is no gauge-invariant 
local operator corresponding to the gluon spin. As a 
consequence, the Altarelli-Parisi evolution equation for 
the 1st moment of Aq{x) and Ag(a;) has been justified 
only on the basis of extrapolation or analytic continu- 
ation of higher anomalous dimensions. Our calculation 
shows that, although nonlocal, Eq.Q-® in fact gives 
the gauge-invariant operator definition of the gluon spin, 
which reproduces the Altarelli-Parisi evolution equation 
even in the Feynman gauge. We can further show that 
the final answer for the 1-loop anomalous dimensions is 
exactly the same also in other covariant gauges contain- 
ing an arbitrary gauge parameter ^. Here, the depen- 
dence on the gauge parameter appears in the graph (a) 
and the left graph of (d), but it is precisely cancelled by 
the ^-dependent terms arising from the quark and gluon 
field-strength renormalization. In this way, the gauge- 
invariance of our gluon spin operator with quantum-loop 
corrections is established. 

In conclusion, we have calculated the 1-loop anomalous 
dimensions of the gluon spin operator appearing in our 
gauge- invariant decomposition of the nucleon spin [l3l. 
We find that the answer is in fact independent of the 
choice of gauge. The answer obtained in the Feynman 
gauge (or any covariant gauges) just reproduces the an- 
swer obtained in the LC gauge, which is also the answer of 
the famous Altarelli-Parisi method. The key factor lead- 
ing to the correct answer is the proper sum over the gluon 
polarization states appearing in the Feynman rule. Note 
that such delicacy does not arise in the Altarelli-Parisi 
method, in which only real splitting processes containing 
physical giuons alone (not virtual gluons) come into the 
calculation. At any rate, we have now confirmed that 
the gluon spin operator appearing in our nucleon spin 
decomposition provides us with the long-coveted opera- 
tor definition of gluon spin, which precisely reproduces 
the same evolution equation as obtained by the Altarelli- 
Parisi method in an arbitrary gauge. This means that, 
in combination with the previous two papers, we have 
solved the long-standing puzzle concerning the gauge- 
invariant decomposition of the nucleon spin. 
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FIG. 1. Momentum space vertices for the gluon spin operator. 
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FIG. 2. One-loop graphs contributing to the anomalous dimensions of quark and gluon spin operators : (a) /^'jqq , 
(c) Ajq^, and (d) Aj^q. Graphs with external self-energies are not shown in the figure. Graphs that are not symmeric with 
respect to the vertical lines through the operator vertex have to be counted twice. 
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